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ABSTRACT 

The cosmological origin of cosmic gamma-ray bursts is tested using the 
method of peak ahgnment for the averaging of time profiles. The test is apphed 
to the basic cosmological model with standard sources, which postulates that 
difference between bright and dim bursts results from different cosmological 
red-shifts of their sources. The average emissivity curve {ACEhright) of the group 
of bright BATSE bursts is approximated by a simple analytical function, which 
takes into account the effect of squeezing of the time pulses with increasing 
energy of photons. This function is used to build the model light curve for 
ACEdim of dim BATSE bursts, which takes into account both the cosmological 
time- stretching of bursts light curves and the red-shifting of photons energies. 
Direct comparison between the model light curve and the ACEdim of dim 
bursts is performed, which is based on the estimated probabilities of differences 
between ACEs of randomly selected groups of bursts. It shows no evidence for 
the predicted cosmological effects. The 3a upper limit of the average red- shift 
Zdim of emitters of dim bursts is estimated to be as small, as ~ 0.1 — 0.5, which 
is not consistent with values ~ 1 predicted by the known cosmological models 
of gamma-ray bursts. 

Subject headings: cosmology: theory-gamma rays: bursts 
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1. Introduction 

The cosmological model of cosmic gamma-ray bursts is commonly accepted, as one of 
the most promising concept of the origin of gamma-ray bursts (GRBs). However, it has 
not been finally approved yet by the observational data. Two critical tests were suggested 
to verify the basic model with standard cosmological sources: dim bursts have to be 
time-stretched and red-shifted in comparison with bright events. 

Gamma-ray bursts are known to have very different time histories, and one hardly 
could check the cosmological effects by direct comparison between particular events. 
These tests should be based on some averaged time-based and spectrum-based signatures, 
which represent the basic properties of GRBs. Several statistical tests have already been 
implemented to compare different groups of GRBs and to resolve the predicted cosmological 
effects. In the case of time-dilation, two scientific groups have checked the average emissivity 
curves (ACEs) derived from the peak alignment averaging of bright and dim bursts, and 
came to opposite conclusions: time-dilation of dim bursts was seen by one group (Norris et 
al. 1994a, 1994b, 1996; Bonnel et al. 1996) and it was not seen by another one (Mitrofanov 
et al, 1992a,b, 1994, 1996). Possible reasons for disagreement were discussed (Band 1994, 
Mitrofanov et al. 1996), and the tentative conclusion has been drawn that the claimed 
dilation of dim bursts was possibly resulted from some systematic in separation of bright 
and dim groups of events. 

On the other hand, in the case of spectral red-shift, all groups involved have a consensus 
that bright GRBs have larger averaged spectral hardness than dim events. It was named, as 
the effect of hardness/intensity correlation (Mitrofanov et al. 1992b, c, 1994, 1996; Paciesas 
et al. 1993, Norris et al. 1994a,b). While this effect was originally seen for the average 
hardness ratio (defined as the ratio of counts at high and low energy channels), recently 
it was found also for the average peak energy < Ep > of vFy energy spectra (Mallozzi et 
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al. 1995). The average peak energies < > of the integral spectra of photons were found 
to correlate with the photons peak fluxes -F^^^^^ at the 256 ms time scale. The effect of 
hardness/intensity correlation may be interpreted, as a result of cosmological red- shift of 
dim gamma-ray bursts in respect to bright events. The corresponding cosmological red-shift 
factor about 1.6-2.2 (Mallozzi et al. 1995) is consistent with original cosmological models 
based on the interpretations of log N - log F distribution (e.g. see Emslie and Horrack 
1994). 

Therefore, there is a discrepancy between different groups about time-dilation of dim 
GRBs in respect to bright bursts, but, on the other hand, there is an commonly accepted 
agreement between them for the hardness/intensity correlation. 

Separate pulses of GRBs are known to squeeze with increasing energy of photons (Norris 
et al. 1986, Fenimore et al. 1995), and, therefore, the average curve of emissivity becomes 
narrower at higher energies (Mitrofanov et al. 1996). According to the cosmological model, 
when bright and dim bursts are detected at the same energy band in the observer frame 
of reference, their time profiles were actually emitted at less hard and more hard energy 
ranges at comoving frames of reference, respectively. Therefore, making a comparison of 
bright and dim bursts from sources with small and large red-shifts, one should suppose that 
the intrinsic squeezing of the light curves of dim bursts due to the increase of energy of the 
emitted photons could partially compensate their stretching due to the cosmological time- 
dilation. 

This paper provides a test of the basic cosmological models of GRBs assuming them 
to be standard sources. It uses the average emissivity curves for groups of bright and dim 
bursts and takes into account the effects of cosmological time- stretching in the observer 
frame together with the internal energy dependent squeezing of bursts light curves in the 
comoving frames. 
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2. Analytic Approximation of the Average Curve of Emissivity for Group of 

Bright GRBs 

The average curve of emissivity (ACE) of GRBs was introduced (Mitrofanov et al. 
1994, 1996), as a general signature of bursts time variability. To build an ACE, all time 
histories of averaging bursts should be normalized by peak numbers of counts Cmax, then 
should be aligned at their peak bins tmax and then should be averaged at each of another 
bins. Comparison between the First, the Second and the Third BATSE Catalogs (Fishman 
et al. 1994, Meegan et al. 1994 and Meegan et al. 1995) has shown that ACE has rather 
stable shape: it has one peak profile with steep rise front and gentle back slope, and its 
width decreases with increasing energy of photons used for averaging (Mitrofanov et al. 
1994, 1996). 

For the present analysis the DISCLA data were used from the large area BATSE 
detectors (LADs) with 1024 ms time resolution on three discriminator channels, number 
1 (25-50 keV), number 2 (50-100 keV) and number 3 (100-300 keV). Two basic intensity 
groups of BATSE GRBs were selected from the Third BATSE Catalog (3B) (Meegan et 
al. 1995): 296 bright bursts with -F^j^af^-* > 1 photons cm~^ s"^ and 332 dim events with 
^max^^ < 1 photons cm~^ s~^. Only bursts with igo > 1.0s were taken into account for 
consideration. 

The group of bright bursts is used as the reference sample to find the analytical 
approximation of ACEhright at different discriminator channels (Figure 1). The function 

''0 I '^max I 

approximates ACE[*^^^^^ profiles at each discriminator channel i = 1, 2, 3 with different 
power indexes aj^p at the rise front (RF) t < tmax and a^^ at the back slope (BS) t > tmax, 
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respectively. Instead of three different functions (1) for each of three channels, a single 
function fbright{t, E) could be implemented, which approximates the shape of ACEi^right at 
different energies E, which correspond to these channels 



f, , (i Tp\ — ( hrightjE) \aRp(E),aBs(E) (r,\ 
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where the functions 



t^.^gHm = • (E/173 keVr (3) 

aRF{E) = afl ■ {E/173 keVY' (4) 
aBs{E) = afs ■ {E/m keVy^ (5) 

represent the change of KCEi,right shape with energy. A difference between three observed 
ACEj,^jg;j( profiles (Figure 1) and the model approximation (2) could be evaluated using the 
function 



^bright -2-^2^ _9/ A^T:.(i,3) N 



{ACE ^;-g^^ - fbrightjt j, Ei)) 

^ . o-iACEtll^ 
where E^ corresponds to mean energies at three discriminator channels i = 1,2,3 and 
tj corresponds to the time bins of ACE curves from j — —19 up to j — +19. Errors of 
observed ACE profiles were estimated from the sample variance. 

The parameters of approximation 4r!g/it — l-SOials s> ^rf — l-31±ai2) ^bs — l-lOiaog^ 



ai = -0.10 ± 0.16, a2 = 0.06 ± 0.09 and = 0.11 ± 0.08 were estimated from the 

bright 



best fitting of all three ACE[^^. profiles at channels i = 1, 2, 3. This fitting leads to 



the minimum aS^™^^ of Exp. (6), which corresponds to rather small value of the Pearson 
criterion: reduced = 0.66 for 108 degrees of freedom. Therefore, one might conclude that 
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Exp. (2) gives a rather good approximation of the observed ACEi,right profiles for the basic 
group of bright bursts at a broad energy range from 25 up to 300 keV. On the other hand, 
rather small values of the reduced points out that the errors of ACEbright was probably 
overestimated by the sample variance algorithm, or there were some correlation between 
them. 

However, the Pearson criterion allows to determine the confidence region for the 
estimated parameters of the fitting function (2). According to Lampton et al. (1976), the 
confidence region for the significance level A could be determined by the 5-dimensional 
contour Scountur in the 6-dimensional parameter space, which is given by the equation 



Scountur — Sj^ight + xlWi i'^) 

where xiW represents the value of distribution for significance A for 6 degrees of 
freedom. Errors ilcr for each of six parameters, as presented above, were estimated from 
the condition that Exp. (6) for Sbright becomes equal to Scountur when the parameter goes 
up and down from the best fitting value, while another five parameters are used as free 
parameters for minimization. Therefore, each of these 12 points could be interpreted as 
±0" deviations from the minimum point along the axes of corresponding parameter inside a 
5-dimensional contour Scountur- 

Actually, these 12 points in the six-dimensional parameter space correspond to 12 
fitting models (2) of the ACE^^^^^^ curves. Were taken all together, they would present the 
±1(7 corridor of analytical approximations around the best fitting model, which leads to 
Sbright — 'S't^gj^t- ^hc boundary curves of this corridor are presented at Figure 1. One might 
see that all these models provide rather good approximation of all three ACE^'^^^^j profiles 
measured at three energy discriminitors. 
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3. Comparison between the average emissivity curves for different groups of 

bursts 

Particular gamma-ray bursts are known to have very different time histories and 
energy spectra. Therefore, ACE curves could be different for particular groups of bursts 
randomly selected from the total data base. Groups with Nj-ep bursts could be defined, as 
representative samples, provided the differences between their ACEs would be comparable 
with the errors from the sample variance for each group. For smaller groups with < A^ep 
a difference between ACE curves could be significantly larger than it would be expected 
from the sample variance. Therefore, the comparison of ACE of different groups has to take 
into account the actual distribution of differences between ACEs profiles due to a random 
choice of contributing bursts. 

Nobody knows how large is the representative sample of time histories of GRBs, but it 
seems from the comparison of ACE curves for IB, 2B and 3B databases that A^ep could be 
about the presently available number of bursts ~ 10^ (Mitrofanov et al. 1997). As it was 
found there, the Pearson criterion provides a rather sensitive test to measure a difference 
between ACEs for any two groups of bursts, namely groups 1 and 11, at any discriminator 
channel i: 

The magnitude S(^i^ij) was used to compare ACEs profiles for randomly selected groups 
of events. It was found that groups with N increasing from ~ 30 up to ~ 300 become 
more and more representative with respect to the full set. In particular, the probability 
distribution P300 of S^^l-u) at discriminator i — 2 was obtained from 10^ random choices 
of two groups with N=303 among the total 3B set of 638 BATSE bursts (Figure 2). This 
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distribution does not depend significantly on the intensity of selected bursts, because the 
main contribution into S(^j_ij) comes from the actual difference of their time histories. 

Thus, Exp. (8) could be used for direct comparison between ACE profiles of groups 
of bright and dim bursts, and the significance of a physical difference S between them 
could be estimated as the probability of obtaining S greater than S^^u) according to the 
distribution -P3oo('S'(7_/7)) provided by the Monte Carlo random choice test (Figure 2). This 
probability distribution will be used below to compare the analytical model based on the 
ACEhright of bright bursts and the actual ACE^im measured for the group of dim events. 

4. Direct cosmological Test Based on the Analytic Approximation of the 

Average Curve of Emissivity 

The simplest test of cosmological model of GRBs could be based on the standard 
candle assumption, which means that everywhere at cosmological distances all sources have 
the same properties in their comoving frame. This basic version of the cosmological model 
assumes that all groups of bursts, provided would be averaged in comoving frames, should 
have the same ACEs. Therefore, any difference between ACEs of bright and dim bursts 
measured in the observer frame should point out on the cosmological effects. 

Let us assume that the emitters of bright and dim bursts have average red shifts z^right 
and Zdim, respectively. If two standard sources at z^right and Zdim emit bright and dim 
bursts with photons energy Eq and variability time scale tq, they would be detected in the 
observer frame at energies E},right = Eq/{1 + Z},right) and Edim = Eq/{1 + Zdim) and with 
variabihty at time scales nright = 7-o(l + Z},right) and Tdim = to(1 + Zdim), respectively. The 
so-called stretching factor could be introduced 
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y{Zbright, ^dim) — (1 + ^dim) / (1 + ^bright) , (9) 

which equals to the ratio of energies of photons Ehright/Edim and/or to the ratio of variabiUty 
time scales TiUm/Tbright at the observer frame of reference, provided they were the same in 
comoving frames. 

To test the basic cosmological model, the analytical approximation fbright{t,E) (Exp. 
(2)) should be transformed into the model function fdim{t, E) according to cosmological 
red-shifting and time-stretching transformations, which has to represent the measured 
KCEdim profiles for the group of dim bursts. According to the assumption of standard 
candles one should postulate 

fdim{t, E) — fbrightiy, E -Y). (10) 

Using the Exp. (2), one might represent Exp. (9), as the following 



fdimit, E; Y) = { rnHy-E),ass(Y-E)^ (H) 

^ ijbright \^ / I ^ '^max \ 

which could be used either as a function of one stretching parameter y, or as a function of 
two red-shifts Zbright and Zdim- 

Figure 3 presents ACE^*^^ profiles for the basic group of 332 dim bursts from the 3B 
database with F^^^^^ < 1 photons cm~^ s~^ observed at three energy discriminators with 
numbers i = 1,2, 3. Expression (11) provides a trial function for the ACE^*^^ profiles with 
the factor y, as a free parameter. To compare the model with observations, the function 
Sdim could be used similar to S^right (6). Table 1 presents the best fitting values Y* for each 
of the three ACE^*^^ profiles fitted separately, and one more value for the joint fit of all 
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three curves together. The errors of Y* correspond to the range of the best fitting values of 
Y for the 12 different models (11) based on the initial model (2) with ±a deviations of its 
six parameters (see Section 2). 

The values of minima 5*^™"^ for the best fitting parameters Y* are rather large, and 
according to the Pearson criterion, the model of equation (11) does not agree with the 
observed ACE^*^^^ profiles for discriminators i — 1,3 and (1 + 2 + 3). Only in the case of 
discriminator i — 2 the model (11) with y*=0.85-0.87 formally agrees with the ACE^^^ 
profile. Moreover, instead of the expected stretching, all the best fitting factors Y* (Table 
1) correspond to squeezing of ACE^*^^ profiles with respect to the analytic approximation 
of ACE^t^j for bright bursts (Exp. (2)). 

However, the classical Pearson criterion based on the x^-distribution could not be 
applied in this case, because it does not take into account the actual distribution of 
differences between ACEs profiles due to a random selection of contributing events. 

A more accurate test of the basic cosmological model is done below, which takes into 
account the probability distribution of 5'(7_//) resulting from the random sampling of 
BATSE bursts (see Section 3). This test has to provide the upper hmits of Zdim for the 
basic cosmological model with standard sources, which could be deduced from the observed 
profiles of ACEbright and ACE^j^. 

According to this model, a group of bursts with fluxes ~ F corresponds to a deflnite 
red-shift ~ z. While in the Euclidean space there is a flux dilution law ~ R~^, which 
establishes the well-known fiux/distance relation for standard sources, the non- Euclidean 
dilution of fluxes from cosmological emitters is influenced by the effects of photon energy 
red-shifting and light curve time- stretching. 

While each burst has a particular energy spectrum, the average spectral distribution 
could be obtained for any selected group of bursts as well as ACEs were obtained for their 



time histories. According to Band et al. (1993), the energy spectra of BATSE bursts 0(-E') 
could be described by the law 



TP E(2+a) 

<j){E) = A ■ ( • e ^f-fc (12) 



100 keV 
if 



if 



where all energies are normahzed by 100 keV. The BATSE database includes the spectral 
data at 2048 ms time scale, which could be used to find the average spectral parameters at 
peak time intervals. For the group of bright BATSE bursts, the average spectral parameters 
at the peaks are < a >= -0.618, < Epeak >= 329 keV and < /3 >= -3.18 (Mitrofanov et 
al. 1997). 

According to the concept of standard sources, one could use the average spectra of 
bright bursts 0(^^'f'^*)(£'), as a standard distribution of photons for all emitters. Therefore, 
one might derive a universal relation between the observed photon fluxes F and red-shifts 
z of corresponding emitters. For two basic groups of 296 bright and 332 dim bursts with 
average peak fluxes < F^i^ht) >^ 6.15±0.35 photons cm^^ and < Fj^i^^ >= 0.53±0.03 
photons cm~^ s^^, respectively, this relation corresponds to the ratio 



I El 

where 



(16) 
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R 



c 



■ [qoZ+{l-qo){l- ^l + 2zqo)] 



(17) 



{l + z)qiHo 



is cosmological distance to a source, Hq is tiie Hubble constant and go represent tlie type of 
cosmological geometry. The geometry of the Universe with critical density is tested below 
with qo — ao — 0.5. The peak flux (photons cm~^ s~^) was calculated in the 50-300kev 
energy range according to 3B Catalog database. 

Using the average values < F^^^*^ > and < F^'^^ > and the average spectral law 
^{bright) (^j^-^ ^ Exp. (16) could be transformed into the relationship between two cosmological 
parameters: an average red-shift Zdim of emitters of dim bursts and an average stretching 
factor Y between dim and bright bursts. Therefore, the Zdim value could be implemented 
into the model function (11) fdimit, E; Zdim), as a free parameter, to check the consistency 
between the basic cosmological model and observed ACE^*^^ curves for dim bursts. 

To do this, one has to put the Zdim value into the model function fdimit, E; Zdim) and 
to calculate the difference (6) between the model and the ACE^^]„ profile at the energy 
discriminator channel i = 2. The estimated value Sdimi^dim) could be corresponded to the 
probability -P3oo('5'(7-7/)) (Figure 2) to find the difference -S'g-//) equal to this value. The 
integrated probability 



could be interpreted, as the probability that cosmological model with Zdim is consistent with 
observed ACEdim profile. Changing Zdim, one might create this way the probability function 
P30o{zdim) (Figure 4). 

To take into account errors in the parameters of the basic analytical model of 
fbright{t, E), the main theoretical model (11) was used together with 12 additional models 




(18) 
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with ±0" deviations from the best fitting parameters (3). They compose the la corridor of 
models around the medium curve which corresponds to the best one (see Figure 4). 

It was found that the probabihty P{zdim) decreases with increasing Zdim becoming as 
small as the level ~ 3 • 10~^ of 3(T fluctuations at Zdim = 0.07 — 0.09 (Figure 3). Therefore, 
one might consider the value ~ 0.1, as the 3(7 upper limit for average red-shift of emitters 
of the basic group of 332 dim bursts with F^^l'^^ < 1 photons cm~^ s~^. 

Two groups of bright and dim bursts arc used for this estimation which been separated 
by the peak flux F^^^^^ — 1 photons cm"^ s~^. In this case one has the largest possible 
number of events in each sample, ~ 300, with the ratio of corresponding average peak fluxes 
of two samples ~ 12. However, one might suspect that selected sample of ~ 300 bright 
bursts might contain a large deal of bursts at cosmological distances, and, as such, a time 
dilation between bright and dim samples could be difficult to resolve. 

Formally speaking, this statement is not correct: in accordance with the basic 
cosmological model, the increase of Zbright value for the bright group results to more and 
more pronounced cosmological stretching of bursts from the dim sample, provided the ratio 
of their average peak fluxes is fixed. Indeed, the Exp. (16) points out that for a given ratio of 
peak fluxes < F^^^^^^^ > I < -^ma^^ > increase of zi^right from the value leads to increase 
of stretching factor Y {zmghu Zdim) ■ Thus, for the ratio < F^liaht) > / < pi^im) ^ q 
one might find Y = 1.2, 1.6 and 1.8 and Zdim = 0.3, 1.1 and 1.7 for Zbright =0.1, 0.3 and 0.5, 
respectively. 

The found 3a upper limit Zdim ~ 0.08 corresponds to the stretching factor Y ~ 1.05 
and Zbright — 0.03. One might conclude that the basic cosmological model with ~ 300 
standard emitters of bright and dim bursts is consistent with the observed ACEj,right and 
ACEdim curves, provided their red-shifts are Zbright < 0.03 and Zdim < 0.1, respectively. 
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However, even taking into account the argument above, one could apply the proposed 
redshifting technique to perform a more conservative comparison between two samples of 
bright and dim bursts, which could be selected by a more stringent criterion based on the 
slope of logN — logF distribution, and which would be truly isolated one from another by 
the sample of intermediate events in between. 

Let us select two samples of 102 brightest bursts with -F^a*^'**^**^ > 4.0 photons 
cm~^ and 100 dimmest events with F^l^'^'^^*^ < 0.41 photons cm~^ s~-^. The 
brightest sample corresponds to the -3/2 part of the logN — logF distribution (see 3B 
catalog, Meegan et al. 1995). There is about ~ 400 bursts with medium peak fluxes 
in between the brightest and the dimmest samples, and the ratio of the average peak 
fluxes < Fl^"^3htest) > / < F^*™™*"^*) >= 49.8 is as large as possible to imply the largest 
cosmological stretching between them. 

For the new sample of the brightest 102 bursts the analytical approximation (2) 
corresponds to the best fitting parameters, which all agree quite well with the estimated 
±1(7 corridor with the basic sample of 296 bright bursts. The best fitting parameters for 
the ACEbrightest curve are 4rLwest = 1-88 s, a^]^ = 1.37, a^^ = 1.26, ai = -0.25, a2 = 0.065 
and ($3 = 0.075. Similarly to (11), these new parameters could be used to build the trial 
function fdimmest{t, E;Y{zbrightest, ^dimmest)) (H) to fit the ACEdimmest curvc of the Sample 
of 100 dimmest bursts. 

The best fitting values of Y** equal 1.01, 0.80 and 0.88 for ACE^l^^^^^ at the three 
energy discriminators with numbers i —1,2 and 3, respectively. The corresponding values 
of reduced are 2.95, 3.20 and 1.90, respectively. Therefore, the best fitting stretching 
factors ¥** between the samples of the dimmest and the brightest bursts do not manifest 
any stretching. These values are similar to the best fitting factors between the basic 
samples of ~ 300 bright and dim bursts, and they all are consistent with the absence of any 
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cosmological stretching. 

However, to find the upper hmit of the stretching factors between the two samples of 
brightest and dimmest events, one has to compare the trial model fdimmest{t, E; Zdimmest) 
(Exp. 11) with the ACEdimmest curve taking into account the sampling statistics of two 
groups. The probability distribution Pioo('S'(7-//)) has to be used for the two sets of 
~ 100 events (see Section 3). According to Mitrofanov at al. (1997), the distribution 
of -Pioo('S'(/-/7)) will have the same shape as the distribution i-*3oo('S'(/-77)) for sets with 
~ 300 events. Therefore, the value of 5'(/_7/) for the 3(T limit will be about the same. 
However, because of smaller statistics, for samples with ~ 100 events the function (8) has 
denominator in ~ 3 times larger than for samples with ~ 300 events, and therefore, the 
difference between two ACEs profiles allowed by 3fT hmit could be in ~ 1.7 times larger. 

Similarly to the basic case of two samples of ~ 300 bursts, the new samples of ~ 100 
brightest and dimmest events are compared by the proposed technique, when for selected 
values of Zdimmest the probability Pimi.^ dimmest) is estimated (see (18)) to get the found 
difference between the model profile fdimmest{t, E; Zdimmest) and the observed ACEdimmest 
curve at the third energy discriminator channel. The corresponding probability function 
Piooiz dimmest) IS presented at the Figure 5. The 3a upper limit of the Zdimmest value is 0.46. 

Thus, when two samples of the brightest and the dimmest bursts with ~ 100 events 
are compared, no significant increase is found for the best fitting stretching factors in 
comparison with the case of two basic samples of ~ 300 bright and dim bursts. At both 
cases one does not see any evidence for stretching effect at all. Using the samphng statistics 
of bursts, the 3a upper limits are estimated of Zdim for ~ 300 dim bursts and of Zdimmest 
for ~ 100 dimmest events, which equal ~ 0.1 and ~ 0.5, respectively. One could suspect 
that the larger upper limit in the second case results from the smaller sampling statistics of 
groups of ~ 100 bursts, and it hardly provides more evidence for cosmological stretching in 
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comparison with the basic case of groups of ~ 300 events. 

However, formally speaking, one has to conclude that the basic cosmological models 
with standard candles are still allowed for gamma-ray bursts provided they correspond 
to the 3a upper limit z^im < 0.1 for the group of dim bursts with F^^'> < 1.0 photons 
cm~^ s~^ or to the upper limit Zdimmest < 0.5 for the group of the dimmest bursts with 
F4*™*) < 0.41 photons cm s . These limits resulted from the different sampling 
statistics of these groups, and further observations of bursts will allow either to decrease 
these limits, or to resolve the time-stretching effect of dim gamma-ray bursts with respect 
to bright ones. 

Two different average photon spectra with power laws a = 1 and a = 2 were used 
for the test of the basic samples also. At the plane Psool^dim) versus Zdim these models 
correspond to upper and lower lines around the main curve, which was found for the average 
energy spectra (Figure 6). Therefore, the shape of the photon energy spectra does not affect 
significantly the upper limit of Zdim- The upper limits of Zdim could be estimated also for 
different parameters of cosmological geometry. Two curves for chance probability Pmoi^dim) 
were derived for two different sets of cosmological parameters (Figure 7): qo — gq — Q.l 
(open Universe) and go = o"o = 1-0 (closed Universe). One might see that these cases of the 
Universe geometry lead to 3(7 upper limits Zdim ~ 0.08 about the same as the case of flat 
expanding Universe {q^ — Uq — Q.l) . 

5. Discussion and Conclusions 

So, the performed comparison of the ACE profiles for groups of bright and dim 
bursts does not allow Zdim larger than ~ 0.1 — 0.5 for the basic cosmological model with 
standard sources. Moreover, the ACE — based limit of red-shift of dim bursts agrees with 
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non-cosmological models of GRBs in the flat Euclidean space. 

There are two well-known estimations of the red-shifts of emitters of GRBs 
according to cosmological models. The flrst one is based on the average parameter 
< ^ l^rnax >= 0.33 ± 0.01 for 3B data base (Meegan et al. 1995). One should expect to 
have < V/Vmax >= 0.50 for homogeneous distribution of standard sources in the Euclidean 
space. On the other hand, the observed parameter < V/Vmax > is consistent with the 
non-Euclidean geometry of expanding Universe. For distant emitters of dim bursts the 
geometry — based upper limit of red-shift was estimated about 0.5-2.0 (Wickramasinghe et 
al 1993). Taking into account the coupling between the spectral shape and the temporal 
proflles of bursts, Fenimore and Bloom (1995) have obtained much larger upper limit ~2-6. 

Another estimation of cosmological limit of the red-shift was based on the effect of 
hardness/intensity correlation of GRBs. The average peak of uF,^ spectra of dim bursts 
was found to be much softer than the average peak of bright bursts (Mallozzi et al. 1995). 
The corresponding ratio between peak energies of dim and bright bursts leads to the 
spectra — based upper limit of red-shift, which was estimated about ~ 1. 

There is an agreement, at least qualitative, between geometry — based and 
spectra — based upper limits of red-shifts of distant emitters of GRBs. These estimations 
result to Zdim ~ 1 or even much larger. On the other hand, the ACE — based upper limit 
of Zdim ~ 0.1 — 0.5 does not agree with either the geometry — based or the spectra — based 
limits. Therefore, the basic model of GRBs with standard cosmological sources is not with 
all available constraints. This is the main conclusion of the present paper. 

Developing a cosmological model of GRBs, one should postulate some kind of 
2;-dependent property (es) of out bursting sources which could ensure the consistency. 
Generally speaking, z-dependence could be attributed to different properties of bursts 
sources, such as outbursts rate density, bursts luminosity, average energy spectra and 
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average light curves. There is a reasonable consistency between geometry — based and 
spectra — based limits of red-shifts of dim bursts emitters. Therefore, one might not suggest 
any intrinsic 2;-dependence either for the outbursts rate density, or for the energy spectra 
of the emitted gamma-rays, because they would both lead to consistent limits of z for 
the model with standard sources. On the other hand, to make the agreement between 
them and the ACE — based limit, one could postulate some sort of intrinsic evolution of 
outbursting sources which leads to intrinsic squeezing of their light curves with increasing 
red-shifts. There are physical conditions in local cosmological space which vary with z: the 
local density of matter, the local temperature of microwave background, etc., but at the 
present time no one knows how much these conditions could actually influence on bursts' 
light curves, if they could at all. Obviously, a priori there is no physical reason to propose 
this kind of evolution, and it could be considered as a pure phenomenological speculation. 

In addition to ACE — based test, the time-dilation tests should also be done with 
another time-based parameters of bright and dim GRBs, such as pulse width, interpulse 
duration, etc. Comparison of distinct time — based signatures for different intensity groups 
of bursts would allow to distinguish the basic effect of cosmological time-stretching and 
energy red- shifting, which should be identical for all time-energy signatures, from another 
effects resulted from 2;-dependent evolution, which should be different for each of temporal 
parameters. The cosmological paradigm of GRBs could be finally approved at these tests, 
and new knowledge would be obtained about intrinsic properties of close and distant GRBs 
sources in the co- moving reference frames. This studies will be done elsewhere. 
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Table 1. Best fitting factors Y* for ACE; 



Energy range, keV Best fit Reduced P(> x^) DOF 



25-50 0.81±0.02 2.3 9.2 • 10"^ 37 

50-100 0.86±0.01 1.28 0.12 37 

100-300 0.82±0.02 2.0 3.0-10-^ 37 

25-300 0.84±0.02 2.0 < 10"^ 113 
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Fig. 1. — At the left line three viewgraphs for ACE^right profiles are presented for three 
discriminator channels number i = 1,2,3. The best fitting approximation profiles (2) are 
shown for each ACE. At the bottom viewgraph [i = 3) the model for i = 1 is shown 
{dash line) to demonstrate the energy dependence of ACE. At the right fine three zooms of 
corresponding ACEs are presented to show the quality of approximations and the boundaries 
of ±(T corridor of models around the best ones. 

Fig. 2. — The probabihty distribution of S^'^^ values divided by number of ACE bins (38) for 
energy discriminator i = 2. It is provided by 10^ random choices of two groups of 303 bursts 
among the total set of 3B database. The value corresponding to 3a standard deviation is 
shown by dotted-dashed line. One sees that sampling statistics allows much larger difference 
between two sets of bursts than it could be expected from the sample variance for each of 
them. 

Fig. 3.— ACE^^ profiles for 332 dim bursts with F^i°24) ^ ^ photons cm s . The best 
fitting models (11) of ACEs at each discriminator channel are shown by sohd lines. 

Fig. 4. — The estimated probability P3oo{zdim) of consistency is shown between the model 
curve (11) based on the ACEbright profiles and the standard cosmological model and the 
observed ACE^j^ profile at the second (i = 2) discriminator channel. The dashed region 
represents the probabilities for ±cr corridor of models around the best one. The dotted- 
dashed line shows the probability level of standard 3a fiuctuations. 

Fig. 5. — The estimated probability Pioo{zdim) of consistency is shown between the model 
curve (11) based on the ACEbrightest profiles and the standard cosmological model and the 
observed ACEdimmest profile at the third (i = 3) discriminator channel. The dotted-dashed 
line shows the probability level of standard 3a fluctuations. 

Fig. 6. — The same probability Psooi^dim) is shown as at the Figure 4 (solid line) together 
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with two another estimations based on the energy spectra with power law: dashed and dotted 
hnes correspond to ck = 1 and a — 2, respectively. 

Fig. 7. — The same probabihty Psooizdim) is shown as at the Figure 4 (sohd hne) together 
with two estimations corresponded to another models of Universe: dashed and dotted lines 
correspond to go — 0.1 and Qq — 1.0, respectively. 
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